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THE FAINTEST X-RAY SOURCES FROM Z = 0- 8^ 
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ABSTRACT 

We use the new 4 Ms exposure of the CDF-S field obtained with the Chandra X-ray satelhte to 
investigate the properties of the faintest X-ray sources over a wide range of redshifts. We use an optimized 
averaging procedure to investigate the weighted mean X-ray fluxes of optically selected sources in the 
CDF-S over the redshift range z = — 8 and down to 0.5 — 2 keV fluxes as low as 5 x I0~^^ erg cm~^ s~^. 
None of the samples of sources at high redshifts (z > 5) show any significant flux, and at z = 6.5 we 
place an upper limit on the X-ray luminosity of 4 x 10'*^ erg in the rest-frame 3.75 — 15 keV band for 
the sample of Bouwens et al. (2006). This is consistent with any X-ray production in the galaxies being 
solely due to star formation. At lower redshifts we find significant weighted mean X-ray fluxes in many 
samples of sources over the redshift range z = — A. We use these to argue that (1) the relation between 
star formation and X-ray production remains invariant over this redshift range, (2) X-ray sources below 
the direct detection threshold in the CDF-S are primarily star-forming, and (3) there is full consistency 
between UV and X-ray estimations of the star formation history. 

Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: evolution — 
galaxies: starburst — galaxies: active 



1. INTRODUCTION 

We would like to be able to determine directly how much 
growth occurs in supermassive black holes in the early uni- 
verse. Unfortunately, however, only the most luminous 
quasars at z > 6 have ever been detected individually (e.g.. 
Fan et al. 2006; Willott et al. 2010; Mortlock et al. 2011), 
and these sources are so massive that their evolutionary 
histories are far from typical. 

We would also like to know the relative contributions of 
the galaxy and active galactic nucleus (AGN) populations 
to the UV and X-ray ionizing radiation in the early uni- 
verse. Estimates from optical data were that AGNs could 
not account for the required UV ionizing flux at z > 3 
(e.g., Bolton et al. 2005; Meiksin 2005). It is now recog- 
nized that complete samples of AGNs — or at least those 
which are not Compton thick — are most easily obtained 
with X-ray observations. Direct searches for high-redshift 
(z = 3.5 — 6.5) AGNs carried out using combined deep 
Chandra X-ray samples and optical imaging data (e.g., 
Barger et al. 2003; Fontanot et al. 2007) also found that 
the ionization from AGNs was insufficient to maintain the 
observed ionization of the intergalactic medium at high 
redshifts. This result was conclusively established at z > 3 
by Cowie et al. (2009) by measuring the contribution of 
AGNs to the ionizing flux as a function of redshift using 
X-ray, optical, and UV observations. 

However, even if the space-density of high-redshift, un- 
obscured AGNs is not sufficiently high to reionize the uni- 
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verse, it does not rule out the possibility that substantial 
accretion onto supermassive black holes is occurring be- 
hind veils of obscuring material. The only way to deter- 
mine this is by probing the faintest X-ray fluxes. Thus, 
the continued deepening of the deepest X-ray images of 
the sky is essential for studying this important issue. The 
latest of these deepenings is the now nearly 4 Ms of expo- 
sure on the Chandra Deep Field-South (CDF-S; Giaconni 
et al. 2002; Alexander et al. 2003, hereafter AOS; Luo et al. 
2008; Xue et al. 2011, hereafter Xll). The 4 Ms exposure 
is deep enough to detect all sources above a luminosity of 

5 X 10*^ erg s~^ in the rest-frame 3.5 — 14 keV band 
at z = 6. Thus, almost all luminous AGNs, even at very 
high redshifts, will be included in the catalog of directly 
detected sources. However, only one source with a red- 
shift greater than z = 5 has been spectroscopically identi- 
fied (Barger et al. 2003) in the Chandra Deep Field-North 
(CDF-N; Brandt et al. 2001; A03) and CDF-S fields. 

Fainter sources are likely to be dominated by star forma- 
tion, low- luminosity AGNs (LLAGNs), or highly obscured 
AGNs. Such sources should have strong breaks in the rest- 
frame UV owing to the Lya forest of the intergalactic gas 
or to the intrinsic Lyman continuum break and should 
be included in optical or near-infrared dropout samples 
identified from the break properties. We can attempt to 
measure the properties of these fainter sources through 
averaging (sometimes referred to as "stacking") analyses 
of optically or near- infrared selected galaxy samples (e.g.. 
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Brandt et al. 2001; Nandra et al. 2002; Reddy & Stcidcl 
2004; Lchmcr ct al. 2005; Worsloy ct al. 2006; Hickox & 
Markevitch 2006, 2007; Treister et al. 2011). This is of 
particular interest for the highest redshift samples, where 
wc may be probing supermassivc black holes at times very 
close to their formation. However, it also allows us to ex- 
amine the X-ray properties of star-forming galaxies over a 
wide range in redshifts (out to beyond z = 4) and to test 
whether the properties of these galaxies change with time. 

Treister et al. (2011; hereafter Til) used the CDF-N 
and CDF-S data to perform sensitive X-ray stacking anal- 
yses on z > 6 galaxy candidates that could not be detected 
individually to look for detectable X-ray signal from the 
population as a whole. The high-redshift samples they 
used {z ^ 6 from Bouwcns et al. 2006, hereafter Bo06, 
and z ~ 7 and z ~ 8 from Bouwens et al. 2011, here- 
after Boll) were selected with the drop-out technique us- 
ing Hubble Space Telescope (HST) ACS and WFC3 data. 
For the z ^ 6 Bo06 sample, Til found significant detec- 
tions (> 5ct) in both the observed soft (0.5 — 2 keV) and 
the observed hard (2 — 8 keV) bands. The stacked hard-to- 
soft X-ray flux ratio that they found was a relatively high 
factor of ~ 9, which they said could only be explained by 
very high levels of obscuration. Moreover, since the ratio 
was a ratio of the stacks, they inferred that there must be 
very few sources with significantly lower levels of obscu- 
ration. This led them to the conclusion that black holes 
grow significantly at early times but are heavily buried in 
gas and dust. 

These results are rather surprising. Galaxies selected 
by the drop-out technique are generally classic starburst 
galaxies; i.e., they are small, blue, and compact with high 
surface brightnesses (e.g., Stanway ct al. 2004). Mean- 
while, the masses of local supermassive black holes are 
generally correlated with the bulge luminosities of the host 
galaxies (e.g., Ferrarese & Merritt 2000; Kormendy & Geb- 
hardt 2001; Gebhardt et al. 2000; Tremaine et al. 2002). 
At z ^ 3, this would be akin to finding AGNs in Lyman 
Break Galaxies (LBGs), which is not very common (e.g., 
Steidel et al. 2004). Moreover, if there were a lot of gas and 
dust present to obscure the central black hole (Til argue 
that the sources must be nearly Compton-thick along most 
directions), then one might expect the dust to also affect 
the properties of the host galaxy, unless it were somehow 
confined to the region around the AGN itself. 

Indeed, Fiore et al. (2011) conducted an independent, 
non-optimized stacking analysis using the Bo06 sample in 
the CDF-S and did not confirm the Til results. More re- 
cently, Willott (2011) provided an independent critique of 
Til that parallels some of the results of the present paper. 

Here we argue that the Til result is a consequence of an 
incorrect background subtraction procedure. In our anal- 
ysis we find a 2a upper limit on the X-ray luminosity at 
z = 6.5 of 4 X 10^^ erg s~^, which is consistent with any 
X-ray production in the galaxies being solely due to star 
formation. 

At lower redshifts we find significant weighted mean X- 
ray fluxes in many samples of sources over the redshift 
range z = — 4. We use these to argue that (1) the rela- 
tion between star formation and X-ray production remains 
invariant over this redshift range, (2) X-ray sources below 
the direct detection threshold are primarily star-forming, 
and (3) there is full consistency between UV and X-ray 



estimations of the star formation history. 

Wc use a standard Hq = 70 km s^^ Mpc"\ = 0.3, 
f2 = 0.7 cosmology throughout. All magnitudes are in the 
AB magnitude system. 

2. X-RAY DATA 

XI 1 presented the X-ray images from the 3.872 Ms 
total exposure of the CDF-S, which covers an area of 
464.5 arcmin^. They provided a main Chandra source cat- 
alog containing 740 X-ray sources detected with WAVDE- 
TECT at a false-positive probability threshold of 10~^ in 
at least one of three X-ray bands: full (0.5 — 8 keV), 
soft (0.5-2 kcV), and hard (2-8 kcV). The catalog 
also satisfies a binomial-probability source-selection cri- 
terion of P < 0.004, which means tliat tlie probability 
that the sources not arc real is less than 0.004. The 
on-axis fiux limits are » 3.2 x 10~^^, 9.1 x 10~^^, and 
5.5 X lO^^'' erg cm~^ s~^ for the full, soft, and hard bands, 
respectively. 

The CDF-S is now almost twice as deep as the CDF-N 
field, which is the only other extremely deep Chandra field. 
For most of our analysis we will therefore concentrate on 
the CDF-S alone, but in some cases we will also include 
similarly selected samples from the CDF-N, where these 
can provide an increase in the sensitivity. For the CDF-N 
we use the images and corresponding source catalog from 
A03, which contains 503 sources selected in the same three 
X-ray bands described above. 

XI 1 shifted the optical positions in deep i?-band 
imaging data — which were already matched to the HST 
ACS Great Observatories Origins Deep Survey-South 
(GOODS-S; Giavalisco et al. 2004) images— by 0'.'175 in 
right ascension and — 0'.'284 in declination in order to 
match them to the radio positions of sources in the Very 
Large Array (VLA) 1.4 GHz radio catalog {5a) of Miller 
et al. (2008) and N. A. Miller et al. (2011), in preparation. 
They then matched the X-ray centroid positions to the 
optical sources detected in the i?-band image to put them 
on a common astrometric frame. In this paper, we moved 
Xll's source positions back to the GOODS-S ACS frame 
by shifting them by — 0'.'175 in right ascension and 0'.'284 
in declination. This allows us to work with the astromet- 
ric positions of the optical sources and spectral data in the 
GOODS-S samples. 

In Figure 1 we show the X-ray image of the CDF-S. 
With small circles we show the positions of all of the X- 
ray sources from XI 1. With larger circles we show the 
positions of the X-ray sources that also have deep near- 
infrared (NIR) coverage from HST WFC3 imaging (Early 
Release Science or ERS, Windhorst et al. 2011; Cosmic 
Assembly Near-Infrared Deep Extragalactic Legacy Sur- 
vey or CANDELS, Grogin et al. 2011 and Koekemoer et 
al. 2011). Given the X-ray positional uncertainties, we 
identified an X-ray source with a NIR counterpart if an 
F160W band counterpart brighter than F160W~ 25 {5a) 
was within 2" of the X-ray position. If more than one such 
NIR counterpart was within the search radius, then we 
identified the X-ray source with the nearest NIR neighbor. 
As pointed out in Xll (see also A. J. Barger et al. 2011, 
in preparation), nearly all of the X-ray sources have NIR 
counterparts. Finally, we mark with a rectangle the part 
of the GOODS-S region that was uniformly covered with 
ACS. The X-ray sensitivity degrades rapidly at large off- 
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Fig. 1. — CDF-S 4 Ms full-band image overlaid with a circle showing an ofF-axis radius of 6' and a rectangle showing the HST ACS uniformly 
covered region of the GOODS-S (Giavalisco ct al. 2004). The small circles show the positions of the X-ray sources in the Xll CDF-S catalog, 
and the larger circles show the X-ray sources with NIR coverage from the ERS (Windhorst et al. 2011) and CANDELS (Grogin et al. 2011; 
Koekemoer et al. 2011) observations using the F125W and F160W filters on HST WFC3. We refer to the area within both the 6' radius 
circle and the uniformly covered region of the GOODS-S as our core region. 



axis radii as the point spread function (PSF) of the X-ray 

image increases, so we restrict our analysis to sources lying 
within an off- axis radius of 6' (very large circle). The over- 
lap of the 6' region with the uniformly covered GOODS-S 
region is 101 arcmin^. This region contains 360 of the 740 
X-ray sources in Xll and 19,294 of the 33,955 optically 
selected sources given in the v2 catalog of the ACS obser- 
vations of the GOODS-S. In this paper we shall refer to 
this overlap region in both fields as our core region. 

3. TARGET SAMPLES OVER THE Z = - 8 REDSHIFT 

RANGE 

We measured the weighted mean X-ray counts s~^ in 
three types of selected galaxy samples. First, we used 
galaxies with secure spectroscopic redshifts. This is a rel- 
atively small sample, and the sources are bright in the rest- 
frame optical. Thus, if more optically-luminous galaxies 
are X-ray brighter (e.g., Lehmer et al. 2005), then this 
sample may be expected to produce higher weighted mean 
X-ray counts s~^. Second, we used a photometric redshift 
sample to probe fainter galaxies. Finally, we used dropout 
selected samples to provide as complete a selection as pos- 
sible in a given redshift range. 

3.1. Secure Spectroscopic Redshift Sample 

For the GOODS-S, we started with the European South- 
ern Observatory (ESO) master redshift catalog compiled 
by A. Rettura in 2004 and subsequently updated by 
Popesso et al. (2009) and Balestra et al. (2010). This cat- 
alog includes redshifts from Cristiani et al. (2000), Groom 



et al. (2001), Bunker et al. (2003), Dickinson et al. (2004), 
Stanway et al. (2004), Strolgcr et al. (2004), Szokoly et al. 
(2004), van der Wei et al. (2004), Dohcrty ct al. (2005), Le 
Fcvre et al. (2005), Mignoh et al. (2005), Ravikumar et al. 

(2007) , and Vanzella et al. (2008). We added to this cat- 
alog 144 new redshifts for previously unidentified sources 
that we measured with DEIMOS in the Fall of 2010, as well 
some additional redshifts obtained by S. Koposov et al. 
(unpublished) that were presented in Taylor et al. (2009). 
In our compilation we only included sources with secure 
redshifts. In the case of the VIMOS-VLT Deep Survey or 
VVDS (Le Fcvre ct al. 2005), we only included sources 
classified as 100% secure. We included a small number of 
redshifts from papers that did not give a quality flag, but 
in these cases the spectral identiflcations appeared robust 
based on the material presented in the papers. 

In our final compilation, wc have secure redshifts for 
1792 galaxies and 158 stars in the fuh GOODS-S field. 
In Figure 2(a) we plot the spectroscopic redshifts versus 
the F850LP magnitudes for the final compilation (black 
small squares). 1158 of the galaxies with secure redshifts 
lie in the core region. Using insecure redshifts would add 
a comparable number of additional redshifts, but many of 
these are clearly problematic based on our inspections of 
the spectra or on our comparisons of the redshifts given in 
the different catalogs, so we prefer to work only with the 
secure redshifts. 

For the GOODS-N, we used the catalogs of Barger et al. 

(2008) updated with as yet unpublished DEIMOS redshifts 
obtained by us and some additional redshifts from Cooper 
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et al. (2011). The GOODS-N sample has secure redshifts 
for 2987 galaxies and 199 stars. The slightly larger sample 
in the GOODS-N as compared to the GOODS-S compen- 
sates for the shallower exposure in the CDF-N, and we use 
both samples in computing the weighted mean fluxes. 
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Fig. 2. — (a) Secure spectroscopic rcdshift vs. F850LP magni- 
tude for all the sources with secure spectroscopic redshifts in the 
full GOODS-S (black small squares; see text for references). X- 
ray sources are marked with a black enclosing square. The lower 
blue histogram (right-hand vertical scale) shows the fraction of 
sources that have not been identified as a function of magnitude, 
(b) Same for the sources with photometric redshifts from Grazian 
et al. (2006). The colored symbols show the b (blue), v (green), and 
i (red) dropouts from Beck06 (plus Bo06 for the i dropout sample) 
that have (a) spectroscopic redshifts or (b) photometric redshifts. 
The more X-ray luminous broad-line AGNs are not selected by the 
dropout techniques. 



3.2. Photometric Redshift Sample 

For our photometric redshift sample, we restricted to the 
deeper CDF-S data and used the GOODS-MUSIC catalog 
of Grazian et al. (2006). (Dahlen et al. 2010 report on an 
alternate photometric redshift computation based on up- 
dated data, but the actual catalog has yet to be released.) 
The Grazian et al. catalog gives photometric redshifts for 
13,820 sources in the fuU GOODS-S field. They reproduce 
the secure spectroscopic redshifts quite well, though there 



are a fair number of outliers: about 4% of sources lie at 
{z~Zp)/{l + z) > 0.2. Much of this seems to be associated 
with the z = 2 — 3 redshift range. At higher redshifts there 
are few problems, presumably because the identifications 
are more secure due to the breaks. In Figure 2(b) we show 
the photometric redshifts versus the F850LP magnitudes 
for the Grazian et al. (2006) sample (black small squares). 
This sample contains 10,071 galaxies in our core region. 

3.3. Dropout Samples 

We used Beckwith et al. (2006, hereafter Beck06)'s list 
of b (1335), V (328), and i (105) dropouts in the GOODS-S 
and Hubble Uhra Deep Field (HUDF). Bo06 compiled an 
alternative list of 522 i dropouts in these fields plus the 
GOODS-N. The Beck06 i dropout sample overlaps with 
the Bo06 sample but contains some sources that are not 
in Bo06 and vice versa. We will consider both of these 
samples separately in our subsequent analysis. 

We show the sources in these dropout selected samples 
in the redshift-magnitude plots of Figure 2 (colored sym- 
bols) . All of the spectroscopically identified sources fall at 
the redshifts expected from the dropout selection. There is 
a small amount of inconsistency between the dropouts and 
the photometric redshifts, as can be seen in Figure 2(b), 
where a number of the dropouts lie at low photometric 
redshifts. The inconsistent fraction is 68 out of 631 ob- 
jects with photometric redshifts in the b dropout sample, 
33 out of 161 in the v dropout sample, and 6 out of 39 in 
the i dropout sample. The median photometric redshifts 
(which we will use for later plots) are z — 3.74 for the b 
dropouts, z — 4.78 for the v dropouts, and — 5.84 for the 
i dropouts (where the latter contains both the Beck06 and 
Bo06 samples). 

Higher redshift {z ~ 7 — 8) sources fall out of the 
GOODS-S F850LP selected sample. Here we use the sam- 
ple of Boll. Although McLure et al. (2011) criticized some 
of the Boll sources (which they consider to be lower red- 
shift interlopers) and gave their own robust list, we use 
the Boll samples here because they are deeper. However, 
we note that we do not detect any significant mean X-ray 
flux in either sample and that the error limits are similar 
regardless of the choice of sample. 

4. OPTIMIZED AVERAGING PROCEDURE 

A fully optimally weighted determination of the mean X- 
ray flux in a sample from the Chandra data would involve 
using variable elliptical apertures that are dynamically ad- 
justed to match the local PSFs and background. However, 
with optimal weightings, most of the signal — even for a set 
of sources that is uniformly distributed over the full Chan- 
dra field — arises at relatively small off-axis radii (less than 
6'), where we may include most of the counts within the 
PSF with a moderately sized circular aperture (e.g., Allen 
et al. 2004). Even at these small off-axis angles the opti- 
mal weighting of the Chandra data involves several issues, 
most particularly which aperture to choose and how far 
out in off-axis angle to include sources (e.g., Lehmer et al. 
2005; Hickox & Markevitch 2006, 2007; Til). 

For faint sources, the noise level is dominated by the 
background. We may write the S/N in an aperture of ra- 
dius r as proportional to f{r, 9)/ {b°'^ xr), where f{r,9) is 
the fraction of the counts from the source contained within 
the aperture at off-axis radius 9 (often referred to as the 
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enclosed energy fraction), and b is the background per unit 
solid angle (see also Til). We show this quantity for the 
2 — 8 keV band in Figure 3 as a function of 6 (black solid: 0' 
or on-axis; red: 1.5'; blue: 3'; green: 5'; black dashed: 7') 
and r, assuming b is roughly constant within the region. 
As Til point out, the optimal r to maximize the S/N in- 
creases with increasing 9, rising from r — 0.4" at 9 = Q' 
(black solid) to approximately r = 5" at 9 = 7' (black 
dashed). However, as can also be seen from Figure 3, the 
dependence on r at larger 9 is very soft, so choosing a 
smaller aperture radius has little effect on the S/N. More- 
over, using a smaller aperture radius has the advantage of 
minimizing contamination by neighboring sources. 

A very small aperture is optimal close to on-axis, but 
the use of such a small aperture places stress on the as- 
trometry (see Til). We have therefore chosen to use a 
constant r — 0.75" aperture within an off-axis angle of 
9 = 3' and a constant r = 1.25" aperture at larger 9. 
This roughly brackets the optimal radius for 9 = 3' (see 
blue curve in Figure 3 for the 2 — 8 keV band), namely 
r = 0.9" in the 0.5 - 2 keV band and r = 1.2" in the 
2 — 8 keV band. The noise in typical samples computed 
with this procedure only differs from that computed with 
a fully optimized aperture as a function of radius by a few 
percent. 

We include sources out to 9 = 6', which, in conjunction 
with the samples lying within the well covered areas of 
the GOODS fields, is our definition of the core region. We 
stress that the results are quite insensitive to the choice 
of the maximum 9, because the high off-axis sources have 
low S/N in an optimally weighted mean. Decreasing the 
off- axis angle at which we include sources to 9 = 5' (4') 
would raise the noise in typical samples by 5% (15%). 
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Fig. 3. — Relative S/N in the 2 — 8 keV band vs. aperture ra- 
dius shown for five values of 9: 0' or on-axis (black solid curve), 
1.5' (red), 3' (blue), 5' (green), and 7' (black dashed). The value is 
normalized to the peak on-axis value. 

With the aperture specified, we can now compute the 
X-ray counts s~^ as C = (S" — B)/{f{r,9) x t), where S 
is the number of counts in the aperture, B{= Tir^b) is the 
number of background counts expected in the same aper- 
ture, and t is the effective exposure time at the position of 
this aperture. C may be negative or positive. We take the 
error on this quantity to be ii^ = B^-^ / {f{r, 9)xt), since, in 



general, the number of background counts in the aperture, 
B, is high enough to justify the Gaussian approximation. 

For a given sample we form the optimally weighted mean 
X-ray counts s""'^, L, and the corresponding error, E^. To 
do this, we first exclude all sources with a known X-ray 
source within 6" or a bright (0.5—8 keV counts above 1000) 
X-ray source within 13" in the existing catalogs (XI 1 for 
the CDF-S and A03 for the CDF-N) to remove any con- 
tamination from these sources. These two values (6" and 
13") correspond to the 90% and 95% enclosed light radii, 
respectively, at an energy of 5 keV and an off-axis radius 
of 6'. This ensures that the contribution from the ex- 
tended wings of the known sources is less than about a 
count in any aperture, even for the hard band. When av- 
eraged over the ensembles, this level of contamination is 
negligible. Our contamination radii are smaller than the 
values adopted by Hickox & Markevitch (2006, 2007) or 
by Til, but tests with random samples show that there 
is little sensitivity to increasing them, and choosing the 
smallest possible contamination radii minimizes the num- 
ber of sources excluded by this step. The adopted proce- 
dure reduces the core area from 101 arcmin^ to 86 arcmin^. 
We form the optimally weighted mean X-ray counts s~^, 
L = Yi{C I E"^) E"^), by summing over the remain- 
ing sources in the sample. The corresponding error is 
El = (l/S(l/i?2))0.5, 

The final issue is how to compute the background counts 
per unit solid angle, b. There are two possible approaches 
to this. One can compute a local background for each 
source using either a surrounding annulus (Til) or a 
wavelet approach (Hickox & Markevitch 2006, 2007), or 
one can compute the average background by randomizing 
the positions of the sources and measuring the values ob- 
tained in the blank- field apertures generated (e.g., Lehmer 
et al. 2005). Here we use both procedures: We use annular 
measurements as our primary background determination, 
and then we use randomized position measurements to test 
that the procedure is not introducing offsets. 

Before we continue, we note that the background we 
are measuring is composed of several elements: an ap- 
proximately uniform instrumental background, a truly dif- 
fuse component containing contributions from hot inter- 
stellar or intergalactic gas, and contributions from indi- 
vidual sources below the direct detection threshold of the 
deep CDF images. While a local estimate of the back- 
ground contributions from the instrument and the diffuse 
background should provide a good average measure, such 
estimates could be more problematic for the unresolved 
sources, which may be clustered. Clustering could, in 
principle, result in an underestimate in determining the 
background subtraction, since the annulus in which the 
background is measured lies further from the source than 
the aperture in which the signal is measured. However, in 
practice, the surface density of X-ray sources at a given 
redshift is small, and the average contribution of these 
sources to the background is also small. (Contributions 
from sources at other redshifts may be viewed as random 
with respect to the sources whose signal is being mea- 
sured.) 

We may use the results of Section 5.2 to estimate that 
the total contribution of all the sources in our photometric 
redshift sample is approximately 2% of the observed back- 
ground in the CDF-S. In smaller redshift intervals, these 
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contributions are correspondingly less and drop rapidly 
with increasing redshift (see Section 6). Thus, the redshift 
interval from z = 1 — 1.1 produces about 0.001 of the back- 
ground, and that at z = 3 — 3.1 produces less than 0.0003 
of the background. These contributions are too small for 
clustering in the galaxy population to perturb the back- 
ground in any significant way. 

In determining 6, we use an annulus of radius 8" to 
22". This choice of inner boundary avoids any contribu- 
tion from the target, and this choice of outer boundary 
provides a large enough number of pixels. In order to elimi- 
nate any contribution to h from bright neighboring sources, 
we exclude 6" regions around known sources (13" regions 
around very bright sources), and we clip pixels with high 
counts. Considerable caution must be taken in perform- 
ing the clip. Because the average number of counts in each 
pixel is small, the counts distribution is Poissonian. Too 
severe of a clip can introduce a downward bias in the es- 
timation of the average background in the annulus, which 
would result in a spurious signal when that background is 
subtracted from the counts in the aperture. This appears 
to be the cause of the z ~ 6 weighted mean signal found 
by Til. 



counts, while the less severe clipping (blue) does not. The 
4 counts clip used in the present paper is the most strin- 
gent clip possible that does not produce a significant off- 
set. However, we note that the results are not sensitive to 
choosing higher values for the clipping. 

We can see this offset even more clearly when we plot the 
two distributions of weighted mean 2 — 8 kcV counts s""'^, 
L. (Red solid for the Til procedure; blue solid for the 
present procedure.) The offset of the red solid histogram 
(8.1 X 10"'' ± 2.2 X lO"*" counts s~^) is almost identical to 
the weighted mean signal found by Til for the Bo06 sam- 
ple (downward pointing arrow). In contrast, the offset of 
the blue sohd histogram (5 x 10~^±2.2 x lO"'' counts s~-^) 
is consistent with the zero weighted mean signal expected 
for random realizations if the background subtraction is 
correct. 

The error in the simulations measured from the disper- 
sion in the realizations is about 30% higher than the formal 
statistical error of 1.6 x 10~^ counts s^^, suggesting that 
systematic effects may be present. Thus, we allow for this 
in assessing the significance of the detections. 
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Fig. 4. — Two distributions of 2 — 8 keV counts s~l, C, from 100 
realizations of randomly located apertures in the GOODS-S core 
region for a sample with a size equal to that of the Bo06 sample of 
2 ~ 6 galaxies. The red dashed curve shows the results when the 
background pixels are clipped above 2 counts (this corresponds to 
the Til procedure), and the blue dashed curve shows the results 
when the background pixels are instead only clipped above 4 counts 
(present procedure). The solid line histograms show the distribu- 
tions of the weighted mean 2 — 8 keV counts s~^, L. The more 
severe clipping procedure of Til (red) introduces an offset in the 
distribution that is almost identical to the value measured in the 
actual Bo06 sample by Til (downward pointing arrow). 

To test both the present procedure and the Til proce- 
dure, we ran 100 realizations of randomly located aper- 
tures in the GOODS-S core region for a sample with a size 
equal to that of the B0O6 sample of z '-^ 6 galaxies. In 
Figure 4 we show the two distributions of the 2 — 8 keV 
counts s~^, C, that we measured. In one case we clipped 
the background pixels above 2 counts (red dashed curve; 
this corresponds to the Til procedure), while in the other 
case we only clipped the background pixels above 4 counts 
(blue dashed curve; present procedure). The more severe 
clipping (red) produces a systematic offset in the measured 



5. RESULTS 
5.1. Secure Spectroscopic Redshift Sample 

In Figure 5 we show the weighted mean counts s~^ in 
both the 0.5 — 2 keV (red squares) and 2 — 8 keV (blue 
diamonds) bands versus redshift for the core samples with 
secure spectroscopic redshifts. In (a) and (b) we show 
the results, respectively, for the GOODS-S and GOODS- 
N fields, and in (c) we show the results for the two fields 
combined. The individual fields show a broadly similar 
pattern. As has been found in previous work (e.g., Lehmer 
et al. 2005), there is a strongly detected signal in the 
0.5 — 2 keV band out to z = 4 (see also Bomans, Zinn 
and Blex 2011). Beyond this redshift, even in the com- 
bined fields, the signal falls below the 3cr threshold, even 
when only the formal statistical error is used. 

The weighted mean 0.5 — 2 keV counts s~^ in the com- 
bined fields are 5.2±1.0x 10~^ counts s~^ for the redshift 
interval z = 3 — 4 and 7.2 ± 3.0 x 10~^ counts s~^ for 
z — 4 — 5. The weighted mean 2 — 8 keV counts s~^ in the 
combined fields are weaker and only detected out to z = 3, 
reflecting Chandra^s poorer sensitivity at these energies. 

Converting the 0.5 - 2 keV (2-8 keV) counts s^^ to 
flux with an average multiplication of 6.8 x 10~^^ erg cm~^ 
(23.9 x 10^^^ erg cm^^), where the conversion factors are 
taken from Til, we find the flux ratio /(0.5-2 keV)//(2- 
8 keV)=0.91 ± 0.15 at z = - 1, 0.95 ± 0.23 at z = 1 - 2, 
and 0.65 ± 0.32 at z = 2 — 3. If the source spectra can 
be represented by power laws, then these flux ratios cor- 
respond to photon indices of 1.93 ± 0.13, 1.96 ± 0.19 and 
1.69±0.28, respectively, showing that the sources are quite 
soft on average. 
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compare with the GOODS-S secure spectroscopic redshift 
sample in the same redshift intervals. In Columns (1) 
and (2), we give, respectively, the minimum and maxi- 
mum of the redshift interval. In Column (3), we give the 
number of sources in the spectroscopic redshift sample. 
In Columns (4) and (5), we give, respectively, the mean 
0.5 — 2 keV and 2 — 8 keV counts with la errors for the 
spectroscopic redshift sample in units of 10"'' counts s~^. 
In Column (6), we give the number of sources in the pho- 
tometric redshift sample. Finally, in Columns (7) and (8), 
we give, respectively, the mean 0.5 — 2 keV and 2 — 8 keV 
counts with la errors for the photometric redshift sample 
in units of 10~^ counts s~^. 

Figure 6 is similar to Figure 5 in that we see a strong sig- 
nal out to z = 4 in the 0.5 — 2 keV band but no detections 
at higher redshifts. For example, for the redshift interval 
z = 4 — 5, the weighted mean 0.5 — 2 keV counts s^^ is 
1.38±0.62 X 10-^ counts s'^ or /(0.5 - 2 keV)=9.4±4.2 x 
10^^^ erg cm^^ s^^. We also again see a weak 2 — 8 keV 
signal with respect to the 0.5 — 2 keV signal; it is only 
detected in the redshift interval z — 1 — 3, where we find 
a flux ratio /(0.5 - 2 keV)//(2 - 8 keV)= 0.67 ± 0.15 
corresponding to a photon index of F = 1.7. 
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Fig. 6. — Weighted mean X-ray counts s~^, L, vs. redshift for 
the sources with photometric redshifts in the GOODS-S core re- 
gion. The red squares denote the 0.5 — 2 keV band, and the blue 
diamonds denote the 2 — 8 keV band. Error bars are ±1ct. Also 
shown are the weighted mean 0.5 - 2 keV counts s"! of the Beck06 
b, V, and i dropout samples (red open squares; these samples are 
from the GOODS-S and HUDF) and the Bo06 i dropout sample 
(red open triangle; this sample is from the GOODS-N, GOODS-S, 
and HUDF). Sources directly detected in the XI 1 or the AOS X-ray 
catalogs are excluded. None of the 2 > 4 samples has any significant 
signal. 



Fig. 5. — (a) Weighted mean X-ray counts s~ , L, vs. redshift for 
all sources with secure spectroscopic redshifts in (a) the GOOD-S 
core sample, (b) the GOODS-N core sample, and (c) the combined 
core samples. The red squares denote the 0.5 — 2 keV band, and 
the blue diamonds denote the 2 — 8 keV band. Error bars are itlcr. 
Sources directly detected in the Xll or the AOS X-ray catalogs are 
excluded. 



5.2. Photometric Redshift Sample 

In Figure 6 we show the weighted mean counts in 
both the 0.5 — 2 keV (red squares) and 2 — 8 keV (blue 
diamonds) bands versus redshift for the sources in the 
GOODS-S core region with photometric redshifts. We 
give the results in tabular form in Table 1, where we 



While the overall pattern with redshift is similar in the 
spectroscopic redshift and photometric redshift samples, 
the weighted mean X-ray counts s""'^ are more than a fac- 
tor of 2 lower in the photometric redshift sample than in 
the spectroscopic redshift sample. As we shall show in 
the discussion section (Section 6), this is a simple con- 
sequence of the fact that the photometric redshift sam- 
ple contains sources down to fainter optical magnitudes 
than the spectroscopic sample. The photometric redshift 
sample, while containing the spectroscopically identified 
sources as a subsample, probes to fainter optical magni- 
tudes, and these optically fainter sources are intrinsically 
fainter in the X-rays, reducing the weighted mean X-ray 
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Table 1 

Mean X-ray Counts Per Second Versus Redshift: GOODS-S 



Spectroscopic Sample Photometric Sample 



znun A^spcc 0.5 - 2 kcV 2 - 8 kcV iVphotz 0.5 - 2 keV 2-8 keV 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0.0 


0.50 


127 


7.64±0.83 


2.59±1.37 


1170 


1.59±0.28 


1.41±0.46 


0.50 


1.0 


350 


9.81±0.50 


2.55±0.81 


2459 


2.29±0.19 


0.09±0.32 


1.0 


1.5 


190 


9.39±0.68 


3.95±1.11 


1500 


3.35±0.25 


1.50±0.41 


1.5 


2.0 


25 


12.2±1.81 


0.09±2.95 


817 


3.76±0.34 


1.76±0.55 


2.0 


2.5 


38 


9.66±1.68 


4.67±2.78 


835 


3.00±0.33 


1.52±0.54 


2.5 


3.0 


32 


8.00±1.82 


4.04±3.04 


697 


2.10±0.36 


0.15±0.58 


3.0 


3.5 


24 


7.82±2.09 


1.19±3.48 


430 


2.46±0.46 


0.21±0.76 


3.5 


4.0 


14 


9.00±2.43 


-5.86±4.02 


253 


1.90±0.61 


-1.03±0.99 


4.0 


4.5 


6 


-1.64±4.57 


-8.19±7.65 


145 


1.00±0.81 


0.07±1.33 


4.5 


5.0 









89 


1.98±0.986 


0.79±1.62 


5.0 


5.5 









46 


0.50±1.33 


.015±2.19 


5.5 


6.0 


5 


1.19±4.07 


11.6±6.79 


31 


2.02±1.70 


0.46±2.81 


6.0 


6.5 









6 


3.64±3.59 


-3.94±5.83 



fluxes. However, because the photometric redshift sample 
contains a larger number of sources, it identifies more of 

the X-ray background light, as measured by the product of 
the number of sources per unit area and the mean counts 
per source, than the secure spectroscopic redshift sample. 
At redshifts between z = and 1, the contribution to the 
X-ray background light by the photometric redshift sam- 
ple is aboiit twice as much as the contribution to the X-ray 
background light by the secure spectroscopic sample, while 
between z = 2 and 3, it is about six times as much. 

The mean counts per source in the photometric red- 
shift samples do not change rapidly as a function of red- 
shift. Thus, the results are not very sensitive to catas- 
trophic photometric redshift errors. If 4% of the sources 
that the photometric redshift sample placed at 2; = — 1 
were really at z = 2 — 3, then the mean counts at 
2 = 2 — 3 would only change from 2.59 x lO"*" counts s~^ 
to 2.55 X 10~^ counts s~^. At z = 3 — 4, the change in the 
mean counts would be even smaller. 

5.3. Dropout Samples 

None of the dropout samples are detected in either the 
0.5 — 2 keV or the 2 — 8 keV bands at any signiflcant level. 
For the Beck06 b, v, and i dropout samples, we find that 
the weighted mean 0.5 — 2 keV counts s~^ are 4.7 ± 3.4 x 
10"^ coimts s-i for z = 3.74, 4.2 ± 7.0 x 10"^ counts s^^ 
for z = 4.78, and 2.5 ±12.3 x 10"^ counts s'^ for z = 5.84. 
We show these with red open squares on Figure 6. 

For the B0O6 i dropout sample, we find that the 
weighted mean 0.5 — 2 keV counts s~^ is — 10.7±5.3 x 10~^ 
counts s~^ (red open triangle in Figure 6), while the 
weighted mean 2 — 8 keV counts s~^ is —13.1 ± 9.7 x 
10~^ counts s~^. In contrast, Til obtained a 2 — 8 keV 
signal of 88±13 x 10~*^ coimts s"-'^, which is highly inconsis- 
tent with the present limits. As we discussed in Section 4, 
the Til detection appears to be the consequence of incor- 
rect background subtraction. 

At higher redshifts {z ~ 7 — 8; Boll sample), we 



find that the weighted mean 0.5 — 2 keV counts s ^ is 
—11.4 ± 5.1 X 10~^ counts s~^, while the weighted mean 
2-8 keV coimts s~'^ is -6.9 ± 9.7 x 10"^ counts s"!. We 
conclude that high-rcdshift sources are generally too faint 
to be detected in X-rays, even in these extremely deep 
stacks. The typical 2a upper limits on the weighted mean 
fluxes are approximately 7 x 10~^^ erg cm"^ s~^. 



6. DISCUSSION 

The X-ray fluxes probed by the optimized averaging 
analysis are, by construction, bounded above by the flux 
detection limits in the direct X-ray catalogs. As we illus- 
trate in Figure 7, even at high redshifts, this places an up- 
per bound on the observed-frame 0.5 — 2 keV luminosities 
of ~ 10*^ erg s~^. We may therefore expect the sources 
in the averaging analysis to be dominated by star-forming 
galaxies with some contribution from LLAGNs (e.g., Bauer 
et al. 2004). The power law photon indices based on the 
broadband X-ray colors (F = 1.7 — 2 for z = 1 — 3; see 
Section 5.1) arc broadly consistent with this interpreta- 
tion. Swartz ct al. (2004) give an average F = 1.7 for the 
ultraluminous X-ray sources (ULXs) that dominate the 
X-ray contributions in strong star- forming galaxies. How- 
ever, the photon indices of ULXs are poorly determined 
above 10 keV, which adds some uncertainty to the com- 
parison at the highest redshifts (z ~ 3) for which we were 
able to measure the weighted mean observed-frame flux 
ratios, /(0.5 - 2 keV)//(2 - 8 keV). The dominance of 
star-forming galaxies at the low X-ray fluxes is also sup- 
ported by fluctuation analyses, which show that the unre- 
solved flux distribution is consistent with an extrapolation 
of the star- forming galaxy log N/log S (Hickox & Marke- 
vitch 2007; Soltan 2011). 
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Fig. 7. — Comparison of the weighted mean 0.5 — 2 keV lumi- 
nosities (units of 10^" erg s~^) of the sources in the spectroscopic 
(red solid squares) and photometric (red open squares) redshift sam- 
ples with the individually detected sources in the CDF-S (blue di- 
amonds) and CDF-N (green triangles). Also shown are the on- 
axis detection limits for the CDF-S (blue curve) and CDF-N (green 
curve) . 
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Dijkstra et al. (2011) predicted the X-ray sky surface 
brightness (i.e., the X-ray background hght) from star- 
forming galaxies using the star formation history of Hop- 
kins & Beacom (2006) and the star formation rate (SFR) 
to X-ray conversion of Mineo et ah (2011). This ahows us 
to make a very simple comparison of the light measured 
from the low-luminosity X-ray sources in the averaging 
analysis with the X-ray prediction made from star forma- 
tion histories measured at other wavelengths. In Figure 8 
we compare the Dijkstra et al. calculation in the 2 — 8 keV 
band (black solid curve; we assume a photon index of 
r = 2) with the X-ray contributions determined from the 
results of the present photometric redshift averaging anal- 
ysis in the GOODS-S core region. In Figure 8(a) we use 
open squares to show for the 0.5 — 2 keV band the quantity 
A'' X S/{Ax (zi — Z2)), where A'' is the number of sources 
used in the averaging analysis, S is the weighted mean flux, 
A is the observed area in square degrees, and (zi — Z2) is 
the redshift interval. We can combine this with the contri- 
butions from the directly detected low-luminosity sources 
in the XI 1 catalog to determine the value for all sources 
with 0.5 — 2 keV luminosities less than 10^^ erg s~^ (solid 
circles). Figure 8(b) shows the same plot for the 2 — 8 keV 
band. 



Fig. 8. — (a) The 0.5 — 2 keV contributions to the X-ray back- 
ground determined from the results of the present photometric red- 
shift averaging analysis in the GOODS-S core region (open squares). 
These contributions combined with the contributions from the di- 
rectly detected low X-ray luminosity sources in the Xll catalog 
give the total contributions to the 0.5 — 2 keV surface brightness 
from sources with X-ray luminosities less than lO*'^ erg s~^ in the 
0.5 — 2 keV band (solid circles), (b) Same for the 2 — 8 keV contribu- 
tions. In both cases, the solid curve shows the Dijkstra et al. (2011) 
calculation in the 2 — 8 keV band of the contribution of star-forming 
galaxies to the X-ray background (for the assumption of a photon 
index of F = 2, this should equal the contribution in the 0.5 — 2 keV 
band), and the dashed curve shows the contribution of LLAGNs 
with X-ray luminosities less than Iff*^ erg s~^, as computed using 
the model of Gilli et al. (2007). 

The agreement between the averaging analysis results 
(open squares) and the Dijkstra et al. (2011) prediction 
are extremely good. As expected from the discussion of 
the photon indices, the 0.5 — 2 keV and 2 — 8 keV sur- 
face brightnesses are almost identical and match to the 
prediction for most of the z = — 3 range over which 
the prediction was made, though the noise levels are much 
higher in the 2 — 8 keV sample. The only exception is 
at low redshifts (2 < 0.5), where the upper bound on the 
flux excludes many of the star-forming galaxies. When the 
results from the directly detected low-luminosity sources 
in the Xll catalog are combined with the averaging anal- 
ysis results (solid circles), the data also then agree with 
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the prediction at these low redshifts. Integrating through 
the combined results, we find that a very large percent- 
age of the contribution to the X-ray background from 
Lx < 10^^ erg star-forming galaxies comes from low- 
redshift objects: 79% at 0.5 - 2 keV and 83% at 2 - 8 keV 
arise from z <^ 2. We also show the contributions from 
LLAGNs with Lx < 10^^ erg s^^ to the 0.5 - 2 keV 
and 2 — 8 keV surface brightnesses, as computed using 
the model of Gilli et al. (2007). These can be significant 
at 2: < 1 but drop rapidly at higher redshifts. 

We may invert this and use the X-ray data to construct 
the star formation history as a function of redshift from 
X-ray samples (e.g., Norman et al. 2004; Lehmer et al. 
2008). We converted the total X-ray surface brightness 
in the 0.5 — 2 keV band of all sources with luminosities 
less than 10^^ erg s~^ in each redshift interval to the 
corresponding comoving X-ray luminosity density in the 
0.5 — 8 keV band, assuming a photon index of F = 2. We 
then computed the SFR per unit comoving volume using 
the relation in Mineo et al. (2011). For the present data, 
we cannot separate the various contributions to the X- 
ray light, so rather than adopting the Mineo et al. (2011) 
best-fitting linear relation obtained from only the resolved 
sample (their Equation 22), we used their linear relation 
obtained from only the unresolved galaxies (given in their 
Section 8.1): 



SFR(Mq yr-i) = 2.7 x 10-4"Lo.5_8 kev(erg s" 



(1) 



Here the SFR is for a Salpeter (1955) initial mass func- 
tion stretching from 0.1 to 100 Mq. The normalization 
is a factor of 1.4 higher than if we had instead used their 
result for only the hard X-ray binaries. We note that this 
calibration does depend on the star formation history. We 
also note that there have been numerous calibrations of the 
SFR with X-ray luminosity, and these vary by up to 40% 
(Grimm et al. 2003; Ranalli et al. 2003; Persic & Rephaeli 
2007; Lehmer et al. 2010; Mineo et al. 2011). Thus, the 
estimated SFR will be uncertain at least at this level. 

We show the star formation history determined from the 
X-ray data in Figure 9, where we compare it with that de- 
rived by Hopkins & Beacom (2006) using multiwavelength 
data. The shape and normalization of the two curves are 
in very good agreement. The normalization agreement 
is probably overly good, since this could be changed by 
any one of a number of factors, such as contamination 
by LLAGNs, cosmic variance, or the uncertainty in the 
X-ray to SFR conversion factor. As discussed above, the 
LLAGN contribution is probably most significant below 
z — I, where it may result in the X-ray estimate of the 
SFR being high. 
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Fig. 9. — The SFR per unit comoving volume calculated from the 
0.5 — 2 keV sources with X-ray luminosities less than 10''^ erg s~^ 
(red squares with la error bars) vs. redshift. The black solid curve 
shows the Hopkins & Beacom (2006) parametric fit to their star 
formation history renormalized to a Salpeter IMF stretching from 
0.1 to 100 M0. 



With these data we can also address the evolution of 
the relation between the X-ray luminosity and the SFR 
as a function of redshift. This relation, Lx = cxx SFR, 
depends on the initial mass function (IMF), the proper- 
ties of binary stars, and the galaxy metallicities, among 
other things. One might expect that the normalization, 
cx, would be higher in low-metallicity galaxies (e.g.. Book- 
binder et al. 1980; Dray 2006; Linden et al. 2010). There 
is some weak evidence for this (Kaaret et al. 2011). In 
turn, this might suggest that cx could be higher at high 
redshifts. Mirabel et al. (2010) invoke this evolution as a 
possible reionization mechanism. 

Dijkstra et al. (2011) tried to constrain such an evolution 
using the contributions to the X-ray backgrounds. Their 
results could only weakly constrain any evolution in cx, 
but this might be due to the fact that the contributions 
to the backgrounds from the higher redshifts are small. 
We can more directly constrain any evolution in cx as a 
function of redshift by comparing how the X-ray fiux de- 
pends on the UV flux as a function of redshift. Since the 
star formation histories, particularly at the highest red- 
shifts, are generally computed from extinction corrected 
UV luminosity densities, this is the most direct test we 
can make. 

For each galaxy in the GOODS-S core region with a pho- 
tometric redshift, we computed the rest-frame UV magni- 
tude at 2500 A by interpolating between the four ACS 
bands in the GOOD-S catalog and then translating this to 
a bolometric fiux, vf^, evaluated at the redshifted wave- 
length. The rest-frame bolometric luminosity at 2500 A is 
vL^ = Andfiyfi,, where di is the luminosity distance. We 
restricted our analysis to z > 1 so that in our lower red- 
shift intervals the above interpolation is valid. At z > 2.6 
where the observed wavelength lies beyond the reddest 
band (F850LP) in the ACS data, we used the F850LP 
magnitude alone to estimate the rest-frame bolometric lu- 
minosity. 
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Fig. 10. — The weighted mean observed-frame 0.5 — 2 keV flux of 
the sources in the GOODS-S core region with photometric redshifts 
vs. the bolometric flux, vf^, evaluated at a rest-frame wavelength 
of 2500 A. If the X-ray source spectra are described by power laws 
with r = 2, then both axes may simply be converted to luminosities 
by multiplying by 4:nd'^ . (a) Correlation over 2 = 1 — 6 for intervals 
of dz = 1. (b) Correlation over 2 = 1 — 3.5 for a finer gridding. 
(c) Same as (a), but retaining in the averaging procedure the di- 
rectly detected low X-ray luminosity (Lo.5— 2 koV < 10**^ erg s"-*^) 
sources in the Xll catalog. In all panels the flux ranges have 
been slightly offset between redshift intervals for clarity. The red 
dashed line shows the on-axis flux limit for a direct detection in the 
0.5 — 2 keV band. The black solid line shows the linear relation 
/o.5-2 keV = 1-57 X 10-3 ^1/(2500 A) given in Equation 2. The 
black dashed line shows the relation that would be obtained using 
the local calibrations of Minco et al. (2011) for the X-ray luminos- 
ity vs. SFR and of Kcnnicutt (1998) for the ncar-UV luminosity vs. 
SFR, assuming F = 2 and omitting any extinction correction. 



Using our optimal averaging procedure, we formed the 
weighted mean X-ray fluxes in unit redshift intervals 
from z — 1 — 6 and in lyfi, intervals stretching from 
gj,g 2m^^ s"^ upwards in half dex intervals. We 
show the results in Figure 10, where we plot the weighted 
mean observed-frame 0.5 — 2 keV flux against vfjy. For 
r = 2 there is no if-correction, and the X-ray luminosity 
in either the rest-frame 2 — 8 keV band or the rest-frame 
0.5 — 2 keV band can be obtained from 4:7Tdf f 0,5^2 kcV- 
Thus, in any redshift interval the plots may be considered 
as luminosity-luminosity plots. 

Figure 10 illustrates a number of interesting points 
with regards to the X-ray fluxes. First, at any given 
redshift, there is a linear relation between the weighted 
mean X-ray flux and z^/^(2500 A). The linear relation 
only holds until we reach the flux level where we be- 
gin to exclude sources that are in the Xll catalogs (red 
dashed horizontal line). This only affects the high flux 
bins (j//^(2500 A)> 5 x 10^^^ erg cm"^ s~^). Retaining 
in the averaging procedure the directly detected low X-ray 
luminosity (Lo.5-2 keV < 10*^ erg s~^) sources in the Xll 
catalog, as in Figure 10(c), extends the relation for the 
low-redshift sources to higher fluxes. 

The underlying reason for the linear relation between 
the X-ray flux and i^/i/(2500 A) is the linear relation be- 
tween the X-ray flux and the SFR. However, since the 
conversion from the UV flux to a total SFR requires an 
extinction correction, it also shows that the average ex- 
tinction correction is similar for all sources regardless of 
the near-UV flux. 

Second, there is little variation in the normalization with 
redshift. For the full data set we find 

/0.5-2 kcV = 1-57 ± 0.07 X 10-3i./^(2500 A) (2) 

(black solid line). However, fits to individual redshift in- 
tervals give almost identical ratios (1.48 ± 0.09 x 10"'^ at 

2 = 1- 

at z - 



-2, 1.62±0.12xl0--^at z = 2-3, 1.89±0.24 x 10" 



4, and 2.25 ± 0.75 x 10-^ at z 



5), 



suggesting that any increase with increasing redshift is 
small. More specifically, the ratio is proportional to 
€{z)cx{z){l + z)'^^^ , where e(z) is the mean extinction 
correction for the UV flux and the (1 + z)'^~^ factor arises 
from the X— correction for a power law spectrum with 
photon index F. For F = 2 and e(z) constant, the rise in 
cx{z) from z = 1.5 to z = 4.5 would be 1.5 ± 0.5, and for 
F — 1.7 and e(z) constant, it would be 1.1 ± 0.4. 

We may also compare the measured relation between 
f 0.5-2 keV and i//iy(2500 A) with that derived using local 
calibrations. For this derivation we use the Mineo et al. 
(2011) calibration of the SFR with X-ray luminosity given 
in Equation 1, assuming a spectral index F =2 to convert 
the 0.5 — 8 keV luminosity to 0.5 — 2 keV luminosity: 



SFR(M„ yr- 



5.4 X 10"''"Lo.5- 



5-2 kcV 



(erg s 



(3) 



and the Kennicutt (1998) calibration of the SFR with the 
near-UV luminosity 



SFR(M„ yr-i) = 1.17x 10"'^3i^L^(2500 A)(erg s" 



(4) 



These are all computed for the same Salpeter (1955) IMF. 
Since we are only interested in the ratio, we only require 
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a consistent choice of IMF. We combine Equations 3 and 
4 to obtain the relation 

io.5-2 kev = 2.2 X 10-^z/L,(2500 A) , (5) 

which translates to 

/0.5-2 keV = 2.2 X 10-^i//,(2500 A)e(^)(l + zf'^ , (6) 

when we allow for the extinction correction, e(z), to con- 
vert from the UV-derived SFR to a total SFR, and when 
we allow for the if— correction. In Figure 10 we show the 
result (black dashed line) obtained from this relation for 
r = 2 when we do not correct for extinction (i.e., €{z) = 1). 
However, to match the observed flux, we need to assume 
that e(z) ^ 5, which is very similar to values estimated 
for the UV extinction corrections at these redshifts such 
as Erb et al. (2006) who obtained a mean correction of 4.5 
for LBGs at z = 2. For F = 1.7, the required extinction 
correction is reduced to e(z) ~ 3. 

We can see from Figure 10 that we can only marginally 
detect the most luminous star- forming galaxies at 2 = 5— 6 
with the present data. The maximum z//jy(2500 A)= 
4 X 10""'^^ erg cm~^ s"-'^ corresponds to /0.5-2 koV = 
6 X 10"^^ erg cm^^ s^^. A substantial detection of the 
z = 5 — 6 sample using an optimal averaging analysis of 
the present sort would require an increase of about a factor 
of 3 in exposure time. 

At higher redshifts {z = 6 — 8) the 2a upper 
limit on i^/j,(2500 A) corresponds to /0.5-2 keV = 7 x 
10~^^ erg cm~^ s~^. This translates to an upper limit 
on the X-ray luminosity at 2: = 6.5 of 4 x 10*^ erg s~^ in 



the rest-frame 3.75 — 15 keV band. The highest redshift 
spectroscopically confirmed source in either the CDF-S or 
the CDF-N is at z = 5.19 (Barger et al. 2003), and there 
are no known sources at higher redshifts with luminosities 
above the threshold luminosities corresponding to the flux 
limits of the X-ray catalogs (see Figure 7). 

Our upper limit on the X-ray luminosity at z = 6.5 
of 4 X 10"^^ erg s~^ is a factor of 20 lower than the lu- 
minosity given in Til, which was based on a supposed 
detection in the less sensitive 2 — 8 keV band. Our detec- 
tion limit is consistent with the X-ray fluxes in the high- 
redshift sources being solely due to star formation, though 
contributions from LLAGNs could also be present. AGNs 
with luminosities above 10*^ erg s~^ thus appear to be 
extremely rare at high redshifts. 
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